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Our peculiar velocity with respect to the CMB rest frame is known to induce a large dipole in 
the CMB. However, the motion of an observer has also the effect of distorting the anisotropies at 
all scales, as shown by Challinor and Van Leeuwen (2002), due to aberration and Doppler effects. 
We propose to measure independently our local motion by using off-diagonal two-point correlation 
functions for high multipoles. We study the observability of the signal for temperature and polar- 
ization anisotropies. We point out that Planck can measure the velocity j3 with an error of about 
30% and the direction with an error of about 20° . This method constitutes a cross-check, which can 
be useful to verify that our CMB dipole is due mainly to our velocity or to disentangle the velocity 
from other possible intrinsic sources. 

Although in this paper we focus on our peculiar velocity, a similar effect would result also from 
other intrinsic vectorial distortion of the CMB which would induce a dipolar lensing. Measuring the 
off-diagonal correlation terms is therefore a test for a preferred direction on the CMB sky. 



I. INTRODUCTION 

Let us consider the motion of an observer (us) with respect to the cosmic microwave background 
(CMB) rest frame, with peculiar velocity (3. The motion of an observer has two effects: aberration, 
which is the apparent deflection of a light bundle due to the motion of the observer, and a Doppler 
effect on the frequency of the photons. The largest effect is due to a Doppler shift, which leads 
(to order f3 = |/3|) to a large dipole in the CMB, even in a perfectly homogeneous sky. Because 
of the motion of the solar system barycenter with respect to the CMB rest frame, it is usually 
assumed that the observed dipole entirely comes from such peculiar velocity effect, leading to the 
measurements /3 = (1.231 ± 0.003) x 10~ 3 , I = 263.99° ± 0.14° and b = 48.26° ± 0.03° in galactic 
coordinates [1-4]. 

However, it seems challenging to disentangle the effect of the local velocity from other possible 
sources, such as an intrinsic dipole or a secondary effect or a global dipolar anisotropy of the 
Universe. In particular, given the fact that the low multipoles appear to have alignments along 
some particular directions [5] , it would be interesting to check whether the dipole itself may contain 
anomalous contributions in addition to the Doppler effect. For these reasons, it would be useful to 
have an independent measurement of our peculiar velocity with respect to the CMB rest frame. 

We propose to do this using the fact that aberration and Doppler effect distort the primordial 
anisotropies and introduce correlations among different multipoles: for the off-diagonal components 
(different ^'s) the correlations are already at first order in /3; for the diagonal ones (same £'s) they 
are only of second order [6] (although [7] claims that using a cut-sky as opposed to full-sky might 
lead to first order corrections on the diagonal terms). Our main aim in this paper is to show that 
the off-diagonal two-point correlation functions can be used by high-resolution experiments, such 
as the Planck satellite, to measure our peculiar velocity, even without using information from the 
dipole. 

In the previous literature only Ref. [8] has shown that peculiar velocities could be measured 
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using the asymmetry in the location of the peaks of the power spectrum between forward and 
backward hemispheres, achieving a possible detection of P at 2 — 3cr for resolution I — 1000 — 1500. 
In [6] the authors have computed the effect on the ae m , in order to show that the contamination 
on cosmological parameter estimation due to our velocity is negligible. Finally [9] has shown how 
to remove both the aberration and Doppler effects from CMB sky maps: by assuming that all 
of the CMB dipole is due to the relative motion between the instrument and the CMB. one can 
use this fiducial value of our peculiar velocity to "de-boost" the data, a procedure that should be 
preferentially carried out in the raw measured time-ordered data. 

Note that an effect similar to the one we discuss here could be generated also by other effects, 
such as a global dipolar anisotropy of the Universe, provided they induce a large dipolar lensing 
effect on the CMB photons. Generally speaking, therefore, the detection of such a correlation can 
be a measure of the existence of a preferred direction in our Universe. 



II. ABERRATION AND DOPPLER EFFECT 



If n is the direction of the incoming light in a reference frame S at rest with respect to the CMB 
frame, the direction n' observed (for the same event) in a reference frame S' which moves with 
a velocity v relative to S can be calculated by applying a Lorentz transformation to the vector 
V = — n, which locally describes the velocity of the light bundle. If we denote by V = — n' the 
velocity in the boosted frame, then the Lorentz transformations we are interested in are (we follow 
a derivation similar to the one in [6]) 

V u - P , Vw 

" 1 - PV U ' Vw 7 (1 - PV U ) ' 1 ] 

where |v| = (3 and j~ 2 = 1 — p 2 . The two vectors u and w are parallel and orthogonal to v 
respectively. They are such that V u = u ■ V = u • n = cos 9 and V w = w • V = — w • n = sin 9, 
where 9 is the angle between v and n. With this notation we can now explicitly write the vector 
n' as follows 

V u ~ P > V w [n - (n ■ u)u] 

1 - /3V U U 7 (1 - PV U ) ^1 - (n • u)2 
= n- u + /3 [n - (n • u)u] 

l + /3n-u 7 (1 + /3n • u) ' 1 ' 

Let us choose now spherical coordinates, using u as the z axis, and 9, ip as the usual spherical 
coordinates for the frame S. We do the same for the frame S", calling 9' the angle between v and 
n'. The observer S' finds a deviation in the direction of arrival of the photons with respect to S. 
The aberration angle is defined as a = n' — n. Its projection along the z axis is given by: 

. /3sin 2 fl 

a ■ z = , (3 

l + /3cos(9 ' 1 ' 

which means that at lowest order in /3 we simply have: 

9-9' = 13 s\n9 + 0{!3 2 ). (4) 

The peculiar motion of an observer with respect to a given source induces also a change in the 
frequency of the observed photons. According to Special Relativity the frequency v' in the boosted 
frame is related to the frequency v measured at rest by the following relation 

v' = 2/7(1 + ph ■ u) . (5) 



This is the usual Doppler effect. 
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III. EFFECT ON THE MULTIPOLES 



It is instructive to consider, first, only the effect of the change of direction - i.e. the aberration 
- on the CMB. This effect can be treated with the same formalism which is used for CMB lensing 
so we will sometimes use the term "lensing" to refer to the effect of aberration. In fact, both effects 
amount to a redefinition of the direction of arrival of the primordial photons. In keeping with the 
standard early universe cosmology, we assume the primordial fluctuations to be Gaussian. It is also 
clear from this treatment that any large scale dipolar gravitational potential would induce a similar 
effect via lensing. 

The observed temperature fluctuation A = AT/T due to lensing is usually computed in a gradient 
expansion [10, 11] as a function of the unlcnscd fluctuation (or primordial, superscript P), where 
the total fluctuation is given by 

A(n') ~ A( p )(n) + A. ( f ) (n)9 ;l (fi) + Ag ) (n)G^ (n)9 ;I (n) + ... , (6) 

where the semicolon stands for the covariant derivative (in the 2-sphere) with respect to the i-th 
coordinate and here V8 = a. Here and in the following we adopt the exact all-sky formalism, as 
opposed to using a flat-sky approximation (which would not be fully adequate here). 1 Note that 
we are allowed to use a gradient expansion only if we look at angular scales larger than the mean 
deviation angle. Since in our case the deviation angle is proportional to /3, this means that our 
treatment is correct for multipoles £ < 1//3 ~ 1000. On smaller scales an exact treatment would 
be required; this is discussed in the case of weak lensing e.g. in [11]. While we do not attempt 
an exact treatment, the perturbative expression is likely to give a meaningful order of magnitude 
estimate, as it happens in the case of lensing in [11]. 

Given the temperature anisotropy A(n) and the lensing profile 0(n), we will need their spherical 
harmonic decompositions, defined respectively as: 



a en 



J dh A(n) y/ m (n) , b tm = J dh 8(n) F/ ro (n) . (7) 



This leads to the standard equation for the lensed ae m : 

ai m - a lm +a lm 

(P) 

where a) is the unlensed signal and the correction is given by 



(8) 



VI" 

where we have used the fact that b^ m is proportional to 5 m o and we kept only the first 2 terms in 
Eq. (6) (and thus Eq. (9) is only correct at order /?). Indeed, 0(n) depends only on the 9 angle 
if we take the axis of decomposition of the ag m in the direction of the velocity v. We have also 
introduced the Gaunt integrals, given in terms of the Wigner 3-j symbols (see e.g. [13]) as follows: 

„ 2 m 3 _ t /(2li + l)(2l 2 + l)(2l 3 + l) ( £ x l 2 4 \f li h la 



u h e 2 i 3 -\J 4?r ^ J \ mi m 2 m 3 J ' {W) 

Moreover, at order (3 the aberration term bi»o just reduces to the dipole, 

& 10 = 24& , (11) 



Although one could think that assuming a flat-sky would be a very good approximation for small scales (say, 
I > 50), the fact that the aberrated power spectrum is given by a convolution of the lensing and un-aberrated 
power spectra propagates the flat-sky errors to high values of i - see [12]. 
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leading to: 



a, 



(-1) (j <fl a e> m b w, (l^j 



4 ro = &io§V^ £+2 ) 





Because of the properties of the Wigner 3-j symbols, we may write this in the following general 
form: 

J L ) ~ r - J p ) + r + JP) (n) 
The Gaunt integrals in Eq. (12), simplify considerably leading to the following coefficients: 

1 v / f 1 

c - = - b ™W^- *>V 4^fr = -^-^V^fr- (14) 

As it has been pointed out by [6], there are also additional terms in this expression due to the 
Doppler effect. We summarize the calculations performed by [6] in Appendix A, which lead to: 

c+ m = p(£ + 2-d) 

^- 1 + rf) V^T' (15) 

where the term due to d contains the pure Doppler effect. Setting the value d = 4 we would 
recover the coefficients obtained by [6] for the multipoles of the intensity integrated over all the 
frequencies. However as discussed in Appendix A, we are interested experimentally in a slightly 
different quantity, namely the Thermodynamic Temperature T. 2 This is derived experimentally 
from the Intensity at some given frequency v via the transformation Eq. (A3) and transforms under 
a boost as in eq. (A4). In this case the above coefficients take the value d = 1, in agreement 
with [14]. From now on all our results will be valid for the Temperature fluctuations. 

Note also that the Doppler term is subdominant for large £; however in an observable quantity 
such as the two point correlation function, the two coefficients c[ m and cj m appear in a combination 
which leads to a partial cancelation in the leading term, proportional to £. In the end, as discussed 
in section IV, it turns out that both effects - aberration and Doppler - contribute in a comparable 
amount. 

The same coefficients can be computed also for E and B polarization fields. We can summarize 
the general result (see again Appendix A) at order j3 

x _ (Px) , -x (Px) , r +x (Px) (K] 

where the index X represents T (temperature), E (electric component of the polarization) or B 
(magnetic component of the polarization). The T coefficients have been displayed above, in Eq. 
(15), while the coefficients for E and B are given by: 



+ E - B Ric + o rn ((l + i) 2 -^ 2 )(l 2 -4) h7 . 



(£ 2 -m 2 )(£ 2 -4) 



2 Note in fact that a subtlety of this analysis is that, despite the fact that / = T 4 , this does not mean that AI/I = 
4AT/T. In fact this is true only at first order, while at second order we need to include AI/I = 4AT/T+6(AT/T) 2 . 
Now, since in the boosted frame T = T (l + j8 cos(0))(l + AT( F )/T), then when taking the square in (AT/T) 2 
we also get a term which goes like (3 cos (e)(AT< p )/T) which represents an additional mixing between neighboring 
multipoles leading to different coefficients c^ m and cj m for the Intensity and the Temperature. 
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where, again, d plays the same role as in the Eqs. 15. For large t they reduce to exactly the same 
coefficients of Eq. (15). Finally it is interesting to note that [6] also found a (weak) correlation 
between the E-mode and the B-mode of the polarization and between the Temperature and the 
B-mode already at 0(f3). However, as discussed in Appendix A, when looking at quantities with 
dimensions of Temperature (instead of Intensities) we find that these correlations exactly disappear: 
namely the Doppler corrections exactly cancels the Aberration. 



IV. MEASURING /3 



As we have seen in the previous section, a mixing between multipoles £ and £ + 1 is present at 
order /3, so that the observed multipoles are corrected as in Eq. (16). We discuss now separately 
the possibility of detecting the aberration signal in temperature and polarization maps. 



A. Temperature 

We can try and construct a non- vanishing two point function which contains this effect. As 
discussed in Sec. I, the most straightforward one, the power spectrum, is unaffected by this since 
it is diagonal in £ space. Although there is actually an effect, it is only at order f3 2 [6], which is too 
small to be observable. 

The off-diagonal part of the two-point correlation function leads to a much larger signal, of order 
/3, instead of f3 2 . We define the following basic quantities: 

Fe m = fl( m %l m • (19) 

These quantities are not rotationally invariant, since they depend on the axis of decomposition of 
the ai m 's. We assume here to know already the direction of the velocity of the observer v of the 
velocity and we choose it to be the axis of decomposition z. We can in this way detect the modulus 
(3 of the velocity. We discuss in section V how to detect the direction of the velocity itself. 
We may analyze the two point function above obtaining (at order /?): 

_ (p)* (p) _ (p)* (p) + (p)* (p) 
rim — a tm a l+lm -\- c e+lm a im a lm + c e, n a e + i m a i + i m 

+ c- a (P) *a {P) + c+ a {P) *a (P) (20) 
Taking now a statistical average we get simply: 

(F lm ) = C- +lm {a {P) n a { /2) + 4m{ a £+*m4+lJ = C i+lm C t + 4m C i+l > ( 21 ) 
( P 1 ) * ( P) 

where Cg = (cif m ag m ) ~ (a) m a \m) ' l& the multipole power spectrum. From Eq. (21) one sees that 
for a static observer (Fi m ) is equal to zero. 3 Note also that this effect does not have an analogous 
counterpart in usual weak lensing. In fact, in lcnsing studies the potentials have vanishing statistical 
average, so that the first nonzero effect arises only at the quadratic level in the potentials. Here 
instead we have used the fact that the velocity has a nonzero known value (j3) = j3 ^ 0. 
For the following it will be useful to work in terms of real quantities: 

n{Flm + Fi-m) = fim: 7T.{Flm ~ Fl-m) = 9lm ■ (22) 

Since fi m — /^_ m , for each £ we have £ + 1 independent real numbers fi m (£ > m > 0), and £ 
independent real numbers ge m - While we have no prediction for the gg m (i.e. their mean value is 
(gim) = 0)j we do have one for the /^ m : 

fim = {fim) = ( c f + l m 8 «m fl <m + i fl <+l m««+l m) = %lm^ + C £m^+1 ' ( 23 ) 



This is true in linear cosmological perturbation theory. However, non linear effects may induce non zero correlations 
between different multipoles. ft would be therefore interesting to compare the correlation patterns emerging in 
this case with the correlations studied in the present paper. 
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In the Gaussian hypothesis, the fg m have at lowest order a diagonal correlation matrix: 

{fimfi'm') — (F im Fti m i) ps -C£Cg + i5a>5 mm >(l + S m0 ) , (24) 
so that the total signal is simply the sum over all the £, m's. Each fi m has a cosmic variance given 

find 



by a £ m = (fL) ~ (fim) 2 , that is 



°L = \c t C l+1 (l + 5 m0 ) - iff*) 2 « \c t C l+1 {l + 5 m0 ) (25) 
(since f£3 is of order f3). We can now build a signal-to-noise ratio summing over all £, m: 4 

Here, and in the following, we adopt the convention that a sum over £, m is to be performed first 
over the range m £ (—£,£) and then over I € (2,£ max ), except where otherwise indicated. If an 
experiment finds a signal with a (S/N) larger than, say, 3, we can safely assume that a detection 
has been made. 

For an ideal experiment, this signal-to-noise ratio, including lensing and Doppler effects simulta- 
neously, is shown by the upper curve in Fig. 1 (both panels). We have fixed j3 = 1.23 x 1CP 3 . On the 
left panel: the dot-dashed curve is obtained by approximating Cg with Ce+i in Eq. (23); the lower 
dotted curve shows how the signal-to-noise changes by considering only the lensing contribution, 
Eq. (14). The right panel shows the comparison between an ideal and a real experiment. Indeed, 
the solid curve includes the noise as foreseen for the Planck experiment, as well as a sky cut. This 
is made by replacing the C/s in the denominator with the following quantities £^ (see e.g. [15]): 

U = -^—{Ci + N e ) , N t = 9 2 a 2 T e e ^ e o^ s 2 > , (27) 

J sky 

where we take for Planck [15, 16] the sensitivity of the best channel: 

a T = 2 x 1CT 6 , e - 7' , f sky = 0.85 . (28) 

It turns out that an experiment which goes up to l max ^ 1000, such as Planck, has the possibility 
of observing the effect with a high Signal-to-Noise ratio. However we should keep in mind that in 
order to use multipoles I > 1/(3 in a prediction an exact calculation is needed, which goes beyond 
the gradient expansion used in this paper, as discussed in sect. III. 

The best fit value for (3 in a given set of observed a^ s , decomposed along the axis of the dipolc, 
is given by the minimization of the x 2 '■ 



obs _ a fTH 
tm I- 1 Jim 



x 2 (/3) = EWt — ' (29) 

£,m 

where we neglected terms of order 1 at denominator and we have extracted the dependence on /?, 
writing f§3 = j3fj^. This leads to the estimator 

fobs fTH \ I I fTH\2 ' 

We can now use Eq. (24) to approximate the variance of j3 as follows 



p 2 p 2 ~ \ h uu+i \s 



(31) 



Starting from a sum over m > we then derive an expression summed over all m's by using the identity fg m 
and Eq. (24). The same for Eq. (29). 
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Figure 1. Signal-to-noise ratio for non diagonal temperature correlation function assuming /3 = 1.23 x 10 . 
Left: ideal case with full computation (solid), approximation as described in the text (dot-dashed), lensing 
contribution only (dotted). Right: ideal case (dotted) and Planck experiment (solid): f s k y = 0.85 and noise 
as described in the text. The horizontal line represents the detection threshold S/N=3. 



B. Polarization 



We can generalize easily the above consideration to the polarization measurements. The observ- 
able quantities are defined as: 

(F X Y] = , X* Y \ _ -Y ( a (Px)*„(P*)\ + r +X (n [Px) * n (PY) ) 
\ r lm I — \ u £m w £+lm/ ~ c e+lm\ a £m a £m I ^ c f m \ u £+lm a £+lml 

_ _-y r XY , +X r XY fon\ 

where X, Y can be either Temperature (T) or polarization (E or B). 
It is useful again to define real quantities: 

1 i p XY , P XY \_ f XY 1 (pXY P XY \ _ XY (oo\ 

~y\ r lm i r £-rn)—J£m ' ^ I r £rn r l-m) = 9£m • \ asi ) 

Again, for each I we have I + 1 real numbers fg m and £ real numbers gi m . Again we have no 
prediction for the g£ m (i.e. their mean value is (gi m ) = 0), while we have a prediction for the fe. m : 

{ffm) — c £+lm ( ^£ XY + C \m^f+\ (34) 

In addition to Cj T which has been already discussed, there are three nonzero primordial correlators: 
Cf E , Cf E and Cf B : 

C( 5i 1 i 2 5 rnirri2 , 

Cf E Si 1 i 2 S rnim2 , (35) 

^£\l2^ra\m2 • 





E \ 




<imi 


a e 2 m 2 / 




E* 


E \ 




limi 


a £2m 2 ) 




B* 


B \ 






%m2' 





CBB c c 



This leads to the following possible observables: fj^i fjmi ffmi ffm anc ^ ffm • Their specific 
shape depends on the form of the primordial correlators. Note also that: 



{ftm ftm) — 7>Cf (1 + S m0 ) . (36) 



We can then define a signal-to-noise ratio for the channel XY as follows: 

TV j rrXXfYY ■ \°'> 
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Figure 2. Signal-to-noise ratio for non diagonal power spectrum for polarization. The dotted curves refer 
to the ideal case (no noise), while the solid curve gives the result for the Planck experimental setup (see 
description in the text). 



We plot in Fig. 2 the signal-to-noise ratio for the polarization channels: ET, TE, and EE. Note 
also that, since TE and ET are practically identical, we plot only their sum, defined as the square 

root of (jf) TE + (jf) ET - In Fig. 2, the dotted curve shows the situation for an ideal experiment, 
while the solid one includes, in addition to the cosmic variance, the noise of the Planck experiment 
(where we assumed that the dominant noise correlations are diagonal in £) and the fraction of 
the sky. For Cf E , the procedure is the same already described for the temperature, but with ar 
replaced by <je = Note that, since the expected BB signal is suppressed with respect to the 

experimental noise, this channel is highly suppressed, leading to a signal-to-noise of at most 10~ 2 
and therefore not detectable with Planck. 

We may finally sum the four off-diagonal temperature and polarization signal-to-noises, in order 
to obtain a total one: 



S I N = JH^I^xy , X,Y = T,E (38) 

V XY 

We show the result in the left plot of Fig. 3. Clearly, this procedure is not exact, since we are 
ignoring the correlations between different components, which could lead to a degradation of the 
signal. One should carry out a Fisher matrix analysis, as we actually do in section VA. It turns 
out that, since the correlation matrix is diagonalised by relatively small rotations (the elements 
of the correlation matrix are given in Eq. (36)), the above estimate for the total signal can be 
considered quite good, also because the total signal-to-noise is essentially dominated by the TT 
contribution (for which the correlations are small). It turns out that, if (3 is close to its fiducial 
value, j3 = 1.23 x 10~ 3 , the Planck experiment should find a S/N in the range between 5 and 10. 

Remember that the S/N ratio is linearly proportional to (3. Clearly, there could be a discrepancy 
between the fiducial value of (3 measured through the dipole, and the value obtained our procedure, 
namely through the non diagonal correlations of high multipolcs. It is then important to understand 
the precision on (3 provided by our method. In the right plot of Fig. 3, we show the relative error 
for (3, i.e. 813/(3, expected assuming that (3 is actually equal to its fiducial value. For other values, 
just remember that 8(3/13 scales as (3~ l . The plot shows that, within our procedure, (3 could be 
measured by Planck with an error ~ 30% (corresponding to a S/N about 4), whereas an ideal 
experiment going up to t — 2000 the error could in principle be reduce to ~ 7% (corresponding to 
a S/N about 15). One should keep in mind though that as discussed in section III, any estimate 
going beyond £ ~ 1000 would require an exact treatment of lensing which goes beyond our approach 
here. 

Finally, a x 2 (/3) can also be constructed for polarization, in an analoguous way as Eq. (29). The 
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2000 



Figure 3. Left: Total signal-to-noise ratio for non diagonal power spectrum. Right: relative error S/3//3. 
We fixed /3 at its fiducial value 1.23 x 10~ 3 . 



final estimator for /3 is obviously (again neglecting the correlation among channels) 



MEE^S^l (eE^) ■ ( 39 ) 

x t, m ^ ^+1 / \ X t,m 1 

while the final expression in the variance equation (31) remains the same with (S/N) given in Eq. 
(38). 



V. MEASURING THE DIRECTION OF THE VELOCITY 



The above signal-to-noise ratio applies if we know from the beginning which is the axis v of 
the velocity and identify it with z. However we would like to detect simultaneously the modulus 
and direction of the velocity. In general, if we want to work in another frame, with z' oriented 
along a different axis than v, we have to apply a rotation to the af m coefficients, using the Wigner 
rotation matrices Lr , (</>, 9, 7) defined in term of three Euler rotation angles <fi, 9 and 7 (see for 
instance [17, 18]). We call 9, (j> the angles that identify the direction of z' in the original frame 
(the one with z). Therefore the parameters in the Wigner matrices can be written simply as 
Tr mm i (cf) ,6,0), where 7, being the angle between the line of nodes and the x-axis, does not affect 
the direction of the velocity and can be set to 0. The rotated coefficients are thus [18] 

af m = E ^mva'afm' • (40) 
-t<m'<l 

In order to build the estimator we define (we drop from now on the bar for the af m coefficients 
in the rotated frame) 

V sXY — n* X n Y (A-\\ 
r lm = a lm a l+lm+s ■ 

where s £ {—1, 0, 1} (although we will restrict only to s £ {0, 1} later). In Appendix B we find that 

(FI*J) = (3 S [ c *-^ +s Cf Y + (-1)^ + X Cf + \] , (42) 

where f3 s = y/4ir/3 Y* a and the coefficients Cp m ' X are defined in Appendix B. A powerful result 
is that each (-F/^ 1 ^) depends only on a single j3 s . This feature naturally implies that the Fisher 
Matrix can be written in diagonal form, as will be shown below. 
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Since it is convenient to work with real quantities, we define the real and imaginary parts of 



psXY 



iXY — 
£m — 



i 



iXY 
£ m 



sXY _ 



2i 



£ m 



F 



Im 



(43) 
(44) 



These quantities in turn satisfy the following (see Appendix B): 



(f£ Y ) = Rc \p.\ [ c |-^ +s cf Y + (-iy Ce 

(gll Y ) = Im %] [4^+. Cr + (-l)'c, 



\s n -s + X n XY 



s„-s + X n XY 



I m 



e+ii i 



(45) 



and 



rsXY rs' hk\ I r 
Jim J I'm' / \J( 



iXY\ 



'em 



' hk\ 



1 



' HI<\ ' 



9lm Y 9l> m^^ — (slm^Vffl'm" ) ~ O^"" ^^ ss ' (^ m0 ~ ^sO^m-m') , 



f sXY a' HK\ 
J I m 9 1' m' I 



f sXY\ 
Jim 



2 
0, 



(46) 



where in the last equation we generalized allowing also for polarization. 

Below we forecast the precision on the parameter determination for a given experiment, such as 
Planck. 



A. Fisher matrix analysis 



We address here the question of the error in the determination of /3 and the direction angles of 
the velocity with a Fisher Matrix analysis. 

Through the Cramer-Rao inequality, the Fisher information Matrix encodes the minimum sta- 
tistical error one is able to recover on the theoretical parameters of the model. This is done under 
the assumption of Gaussian-distributed signal and noise and that the data sets are independent . If 
£(x|p) is the probability of observing a set of data x given the true parameters p, i.e. the fiducial 
model, the Fisher Matrix is defined as: 

„ / d 2 In C \ , , „. 

\ dpiPj I 



where the average is taken over the set of data. Assuming that the likelihood function C is Gaussian 
and that our observables are well defined by (45), our Fisher Matrix is the sum of two matrices 
Fj 7 - — J~ij + Qij that read: 



iX,Y 



dp* 



/XY 



d Pj 



0* = EE 

I,ra,sX,Y 



(48) 
(49) 



where X = XiX 2 = TT, EE,TE, ET (and so for Y) denotes temperature, E-channel polarization 
and temperature-polarization cross-correlations (neglecting the B-channel polarization set of data) , 
driven by the (4x4)- covariance matrices S g jxy- Notice that in general (f^ E ) ^= (/| r if T )- 
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— TT,TT -TT.EE ^TT,TE -TT.ET 



S t = 



~E] 

"/ 

x 

/ 









i — iEE , EE 


— EE,TE 


"EE.ET 
"} 


~TE,EE 


"TE.TE 

"f 


"TE,ET 
a f 


"ET,EE 


— EE.TE 


"EE.ET 



s n = 



/ ~TT,TT ~TT,EE ^TT,TE ~TT,ET \ 

'—'a 



"EE,TT "EE, EE "EE.TE "EE.ET 

"9 "9 "9 "9 



7 TE,TT 
J 9 



7TE.EE 

J 9 



7TE,TE 

J 9 



7TE.ET 

J 9 



i -ET.TT ^ET.EE ^EE,TE ^EE.ET 
\ "9 "9 "9 "9 



(50) 



where 3 f ' — {fi* lX2 ftm Y2 ) an( ^ ^g"' Y = (dim 1 * 2 9em~ Y2 ) are tne covariance matrix elements as 



written in (46), with noise terms 

eJ T = {cJ T + N[ T )/f sky 



£f E = (cr+Nn/f sky , 



iEE 



EE\ 



(51) 



for temperature, E-channel polarization and their cross-correlation, respectively. The overall un- 
certainty encloses cosmic variance and experimental noise, as in (27), including a correction due to 
the galactic cut. 

The inverse of the total Fisher matrix will eventually be the parameter covariance matrix, whose 
diagonal elements represents the marginalized variance errors on p. 

In this Section it is more convenient to build estimators using cartesian coordinates, in which 
f3 is decomposed into three components {f3 x , f3 y , f3 z } (see Appendix B for definitions and detailed 
calculations). As explained in Appendix B we shall count independent quantities in the following 
way: s G {0, +1} with 



m €{—£,£} for s = 1; m£{0,£} for s — (s = — 1 discarded entirely). 
With this convention the key quantity that enters the Fisher Matrix are 



(52) 



d{9t* Y ) 



\Jtm I >r ~ 

(9t* Y ) r = 



h sXY 



Sal 

V2 



(53) 



■ sXY 



where the coefficients hg* Y are given in (B19) of Appendix B. 
A x 2 can be defined as usual, using eq.(B20): 



xHM = E E 

X,Y i,m,s 



If, 



sXY, OBS 



fir \//m ) r 



sXY, OBS 



£f Y £f Y (l + S s0 S m0 )/2 



E 

£ : m,s 



\9c. 



gXYtfXY 



£fJi/2 



(54) 



where // 7 ? Y ' OB and #f™ r '° B are the observed ones. The Fisher Matrix can then be written as 



EE (fXYcfXY 



/fsXYX /fsXYX 

\Jlm l. r \Hm ly lsXY\ I n sXY\ 

l + o s0 o m0 



EE 



X,Y 



fXYpXY 



, (h1* Y Y 

,m=0 



E 



(S r 



Sry) 



IXY\ 



(55) 



which again is a diagonal matrix. Now, just as we did in deriving (26) we shall extend the first sum 
also to negative m values, which in practice incurs a factor of 2 and the disappearance of the <5 m o 
term. We thus have 



EE 



X,Y 



1 



XYifXY 
-l+l 



W.m Y ) 



+ (Srx + S ry ) (h 



IXY\ 



(56) 
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Finally, from (B19) one can check that: 



E 



1XY\ 
£m J 



E 



(57) 



This implies that F 2 



F zz . In other words, the estimates of the errors along each direction 



are the same and we only need to obtain one single standard deviation. This result could have been 
anticipated on physical grounds, as no direction is preferred in this analysis. 

Consequently, we can also define the estimator for the three cartesian components 



E( F_1 ) 

r' 

EE 

X,Yt,tn,s * £ +! 



,v r , 

2 



F 



(58) 



sXY, OBS 
Im 



I fsXY\ 



EEE 



pXYpXY 

771= — £ I ^+1 



X.Y 



1 + 5 s o8 m o 



Urn 



sXY, OBS / n sXY\ 
km \ytm /,, 



ixy, obs^ , „ijry, obSj 



h 0XY rOXY, OBS r 
"*7i Jem ° r 



(59) 



where in the last line we used (53) and the fact that h®* Y = h®_ Y to extend the s = sum to 
negative values of m. The estimator of the magnitude of the velocity is then exactly 



(60) 



and we recover the previous result if we align the z-axis with the correct direction. In conclusion, 
in cartesian coordinates the Fisher Matrix is diagonal, the predicted limits are independent on the 
fiducial magnitude and direction of (3 and the statistical errors in parameter space are the same 
for each direction (a x = a y = a z = a). From this analysis the angles 9 and 4> can easily be 
reconstructed from cartesian components. 

Hence, returning to the case where the fiducial velocity is aligned with the z-axis, we can relate 
for small angles the error on the absolute value of the velocity with the error in the magnitude of 
the direction: 



(61) 



The final statistical accuracy on 9 and /3 as a function of £ max is plotted in fig. 4 for the TT channel 
and for the polarization channels, TE+ET and EE. The errors in the different channels may be 
added in quadrature, as a first approximation, as in eq. (38). As it can be seen from the figure, the 
direction of the velocity can be measured with a precision of about 20° for Planck using multipoles 
up to I ~ 1000, in good part due to the TT-channel. An ideal experiment with no instrumental 
noise on the polarization channels could instead go down to a precision of about 4° if it goes up 
to I = 2000, although we remind the reader again that any estimate going beyond I ~ 1000 would 
require an exact treatment of lensing beyond first order. 



VI. CONCLUSIONS 



In this paper we have shown that a measurement of our peculiar velocity can be achieved using 
the information contained in the Cosmic Microwave Background at high multipoles £, since the 
velocity of the observer distorts the primordial signal through aberration and Doppler effects. We 
have proposed to measure this distortion as a non-zero signal in a two-point correlation function 
flfai%ira 2 , which couples neighbouring multipoles ^'s. We have found that the correlation function 
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Figure 4. Forecast of the error SO — 5/3/fi (plotted in degrees) in the determination of the direction 
as a function of the multipole I using the TT, TE+ET and EE data, as well as the final error obtained 
(adding in quadrature) all of these both for an ideal experiment (dotted line) and for Planck (solid line). 
No assumption is needed on the fiducial value of /?. Compare with Figure 3. 



diagonal in m space can be used to determine the magnitude of the velocity j3, while the off-diagonal 
one can be used to measure the direction of the velocity and we have forecasted the error on the 
determination of such quantities using the Temperature and Polarization channels both for an ideal 
experiment and for the Planck satellite. 

Going up to £ r-i 1000, the error on the amplitude of the velocity f3 which can be obtained is of 
about 30% — 7% (corresponding to a S/N about 4 — 15) for Planck and for an ideal experiment, 
respectively. The direction of the velocity can also be measured with a precision of about 20° — 4° 
respectively for the two cases. 

The method employed in the present paper relies on a gradient expansion, which is valid for 
multipoles £ < 1/(3 ~ 1000, while it would require a more refined treatment which does not rely on 
a gradient expansion for I > 1000. The precision which can be obtained going up to £ = 2000 can 
be estimated with our method and it leads to an improvement of about a factor of 2 in the ideal 
case, and less than 2 for Planck. 

Our proposal constitutes an important independent measurement of our peculiar velocity, which 
can be compared with the usual measurement, through the Doppler effect on the dipole I = 1. This 
can be used to distinguish peculiar velocity from an intrinsic CMB dipole, provided the latter does 
not induce at the same time a dipolar lensing with the same structure in the off-diagonal correlation 
function. 
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NOTE ADDED 

When this paper was concluded, we became aware that A. Kosowsky and T. Kahniashvili in [14] 
independently obtained results, which are consistent with ours, on the possibility of measuring 
the magnitude of the velocity /3, analyzing the Temperature coefficients. In the revised version of 
the present paper we have decided also to analyze the correlation coefficients for the Temperature 
fluctuations, as in [14], instead of the Intensity, as in [6]. In fact a subtlety of our analysis is that 
the multipoles for the total integrated Intensity and for the Temperature do not just differ for an 
overall constant but they have a difference dependence on f3, due to the different Doppler effect. We 
have updated all the figures in the present version of the paper, using Temperature fluctuations, as 
stressed in sect. Ill and Appendix A, and the results now exactly coincide with [14]. 

Then, while [14] has analyzed only the detection of the magnitude of the velocity using the TT 
correlations, we have also addressed here the issue of the detection of the direction of the velocity 
and we have used also the polarization channels. 



Appendix A: Corrections to af m due to Doppler and aberration 

We now calculate the modifications induced by the aberration and the Doppler effect on the 
quantities aj m , af m , af m , respectively the coefficients of the expansion in spherical harmonics of the 
temperature and of the components of the polarization tensor. To separate the Doppler contribution 
from the one associated to the aberration, we introduce a parameter d, such that d = indicates 
that the Doppler effect has been neglected and d = 1, instead, corresponds to the case in which the 
interplay between Doppler and aberration is fully taken in to account. 



1. Temperature 

Under a frame transformation the CMB brightness transforms according to 

(/ \ 3d 
V -J n) . (Al) 

Now, in current CMB experiments, brightness measurements (including those related to polariza- 
tion) are conventionally translated into an equivalent value for the thermodynamic temperature T 
through the relation (see for instance [19-21]) 

V* X V* 

where the function f{x) has been defined in the last step and a; is a frame invariant quantity given 
by x = v/T. For small frequencies, i.e. x — > 0, the function f(x) — > 1 and Eq. (A2) provides 
the antenna temperature which by definition is equal to the ratio I{v)/v 2 . In general, by solving 
Eq. (A2) one finds 

T= , - - N , (A3) 



which is the temperature actually measured in a CMB experiment. Since the ratio i/ 3 /I(v) is frame 
invariant (see Eq. (Al) with d = 1), from Eq. (A3) one can easily read the transformation properties 
of the temperature under a Lorentz boost, namely 

T'(n')= (^)V(n), (A4) 



where again the parameter d has been introduced to separate the Doppler contribution from the 
one due to aberration. 
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We then start our calculation from the following observation 
T'(n') = $>^,y,, m ,(n') 

I'm' 

id 

Cm 



= [ 7 (1 + n ■ f3)f ^aJ m Y tm (n) 

I III 

= E E j / d ^ M 1 + ' /3)]^; m '(n)^ m (n(n)) U W (A') , (A5) 

where we expanded the function of n', [7(1 + n(n') • f3)] d Y^ rn (n(n')), in the basis Y^' m '(n'). 
This expression leads to the identity 

4L' = E/ [7(1 + A(n) • /3)] d r; m ,(n)r, m (n(n)) 

= E/ d ^/ ^^(A'-n)aL[7(l+n(fi)-^)] d ^w(n)y« m (A(n)) 
= 2/ dn'«L[7(l + fi(fi')^)f^(fi')^(fi(fi')) 

Ira 

= W dhaj m [7(1 + n • /3)]^ 2 r; ro ,(n')F, m (n) 

= J2 I dha L' [7(1 + n ' /3)] d_2 r; m ,(n')^ m '(A) , (A6) 

£ J 

where we used the relation dh — [7(1 + n • /3)] 2 dh' and in the last line we assumed that the relative 
velocity (3 points in the direction along which we expanded in spherical harmonics. In this way a 
Kronecker delta removes the summation over m. 

Since we already know that the aberration is a correction which depends from /3, we expand 
Eq.(A6) in powers of j3 and keep only the leading terms in this expansion. We therefore need the 
following Taylor series (/i = n • f3 / j3 and p! = n' • (3/f3) 



[7(1 + n • P)] d - 2 = (1 - (1 + rf) d 2 * 1 + (d - 2)11/3 + 0((3 2 



Y em (ti) = Y im (ti)\t3=o + ^>=o P + 0((3 2 ) 

= Y im (h) + ^(1 - p 2 )fi + 0{f5 2 ) . (A7) 

By means of the following property of the spherical harmonics 

d 

_ ^~d^ lm = ^(e+i)m¥(t+i)m - (£ + l)H gm Y^_ 1 ) m , (A8) 

where 



H '- = fj^j- < A9 > 

we can now write Yj m (n') as follows 

Y tm (p!) = Y tm (h) - pi m {i+l)m Y (l+1)m {h) -(1+ l)H lm Y { i_ x)m {h) \ + 0{i3 2 ) . (A10) 



1G 



Replacing the previous power series expansions in Eq.(A6) one finds (we renamed the indices 
and £') 

0^(1/) ~ ^ (dn [1 + (d- 2>/3] aJ, m Y e , m (n) 

[y; m (n) - /?(w ( , +1)m y ( * +1)m (n) - (£ + l)If, ro Y ( *_ 1)m (n 



EM 



1 + (d- 2) M /3]a? m y rm (n)y; m (n) 



/3a^ ro (^)^'m(n) 



£iJ (£+1)m y ( * £+1)m (n) - + l)iJ fo „r ( *_ 1)m (n)] } . (All) 
The integral over dn can be done by using the following property of the spherical harmonics 
i/iy, m (n)y; m (n) = J dn|c ( ^ +1)m y ( ^ +1)m (n)F^(n) + Q/ m %_ 1)ro (n)F/ m (n)| 

= Hg m 5{ii + i)i + i?(£+l) m <5(£'-l)( • 

This leads to the expression 



/ 



dn, 



(A12) 



a fm — a Jm ~ /3-ff(f+l)m 



(3H (m -{l + d-l) 



\e-d + 2)]af i+1)m (u') 

0(/3 2 ) 



T 

"(£-l)m 



2 Jm + { c lm a J-lm + c i"m a J+l m) 



(A13) 



2. Polarization 



Analogously to the case of the brightness, also the CMB polarization measurements are con- 
ventionally translated into units of temperature. As a consequence, it is convenient to introduce 
a polarization tensor P ab (h) with components related to the Stokes parameters and transforming 
under a Lorentz boost according to Eq. (A4) . Its spherical harmonics expansion reads as follows [6] 



P°*(n) = - V ( ^ m ~ ia f»J-2^™( n ) 







(af m + iaf m ) +2 Y em (h) 



{(J 2 — l)^s^m — 4ff(i+l)m%l)m ~ {(■ + l)s-ff«m^-l)m + 1(1 -f 1) s ^ ri 



i{t + l) 



(A14) 



were ±2Yi m are the spin- weighted spherical harmonics [22]. 

Expanding in the basis ±-iYi m similarly to what we did for the temperature, one finds 

[af m , ± iaf/ m ,] = Wdn [af m , ± ia? m ,] [ 7 (1 + n ■ 0)f- 2 ±2 Y e *, m ,(ii') ±2 Y em ,(h) . (A15) 
Then, by means of the following properties of the spherical harmonics [6] 



(A16) 



where 
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with s = ±2, and keeping just the leading terms in /3, one finally finds 



and 



E \&t'i - P +2H(£ +1 ) m (£ -d + 2)^(^/_i) - +2 H lm (l + d- l)5m> +1 ) } af,, 
t> 

-iPj2 { 4f^^af, m + 0(P^ 



£{£+!) 



— a Fm 



im u i-lm + c tm a i+lm) (A18) 

^ \&t'l - P +2H(£ +1 ) m (£ — d + 2)^(^_!) - +2 Hg m (l + d— l)(5 f (£/ +1) | af, m 
t> 

nS -^ (2d-2)m 



— a fm 



£{£+!) 
( C em a f-l7 



+ ct B a B 



"lm u t+lm) > (A-19) 

where one sees that for the cross E, B terms the inclusion of the Doppler effect (d = 1) exactly 
cancels that of aberration, resulting in no correlation between both polarization modes. This 
should be contrasted with the results in [6], that found a nonzero correlation between E and B 
modes already at O(p). Here it exactly vanishes because we are considering quantities which scale 
under a Lorentz boost according to Eq. (A4) instead of according to (Al). 



Appendix B: Analysis in cartesian components 

If we take a generic velocity, not necessarily oriented towards the z axis, considering the bi Tl 
coefficients of eq.(7), obtained at lowest order in /?, and rotating them as in eq.(40) we get 



47T 



10 



Y?MP 



-Ps , 



(Bl) 



where s — { — 1, 0, +1} and where in the last equality we defined the new variables {/?-, Po, P+} for 
convenience. Note that from the definitions of p s one also has that the cartesian components of P 
are just 



P x =-V2Re[P+], 
P y = V2Im[P+], 
Pz = Po- 



Now, since 



we can write 



£(-l) m G- 



m m—s s 
I £' 1 



£'(£' + l)-£(£ + l) + 2 ( P ) 



l's 



s=-l 



c s ~a {P) +c s+ a (P) 



with the new coefficients given by 



= (-!)'(* + 2)4 



' (£ + 1 + m - s)(£ + 1 - m + s) (£ - sm + s 2 ) 
{A{£ + l) 2 - 1) (1 + s 2 ) (£ + sm) 



(£ + m- s)(£ - m + s) (I + sm + 1 - s 2 ) 
(4£ 2 -l)(l + s 2 )(^-sm + l) 



(B2) 



(B3) 



(B4) 



(B5) 
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Note that these coefficients exhibit the following interesting property, which can be easily verified 
by inspection: 

4i = ciLi- (B6) 
In order to build the estimator let us define again (cf. (41)) 

Ftm = a *tm a l+lm+s ■ (B7) 

which satisfies F% = -F^_ m . 

Note that we have rotated eq.(13) with the coefficients (14) getting eq.(B4) with coefficients (B5) 
for the aberration effect, finding that s can be only 0, +1,— 1. Including also the Doppler effect 
only amounts to rotate with the coefficients (15), which only differ by a d~ dependent prefactor, so 
that we get a similar result with the new coefficients: 



4i = (-ini + 2-d), 

4-=-(e-i + d)J { - 



(£+l + m- s)(£ + 1 - m + s) [£ - sm + s 2 ) 
J (4(£ + l) 2 - 1) (1 + s 2 ) (I + sm) 

+ m — s)(£ — m + s) (£ + sm + 1 — s 2 ) 



(B8) 



(U 2 - 1) (1 + s 2 ) (£-sm + l) 

Now by applying (B4) to (B7), we find that 

(FU = ft [cfo m+s C t + (-1)^+ C i+1 ] , (B9) 

in which we used that (/3 S )* = (—l) s /3~ s and the usual relation {a*i m a>i'm') — CgSwS mm i . This is a 
powerful result since each (F£ m ) depends only on a single ft. This feature naturally implies that 
the Fisher Matrix is a diagonal matrix. For the 4-point correlation functions, assuming Gaussianity, 
we find that 

(Fl m (Fi: m ,y) - (F/ m )((F/, m ,)*) = C e C e+1 5 u ,S mm ,6 ss , . (BIO) 

Since it is convenient to work with real quantities, we define the real and imaginary parts of F£ m 
by: 



ftm — 2 \Ftm + (Fg m ) ] , 



(Bll) 



These quantities in turn satisfy the following: 



(fU = Re [ft] [4~ lm+s C t + (-l)'cJZ+ C t+ 
(g s im ) = Im [ft] [4+ lm+s C t + (-l) s c7 r ;+ C t+ 



(B12) 



and 



fimft'm'} — {fim^{fl'm>) — ^CtCg+iSw (8 ss '8 mm i + 8 s - s /8 m - m >) , 

[g£ m 9i> m ') ~ (9tm)(9i'm') = \ciCi +1 8 u > {6 ss i5 mm i -8 s - s ,8 m - ml ) , (B13) 
^//mSl'm') — W " m^\£ 'I' m') = 0- 
A crucial property regarding these quantities is 

// m = (-l)'/r-m and 9t m = (-l) s+1 9il m - (B14) 
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This implies that all jf m and g\ m for negative m values can be written in terms of their counter- 
parts with positive m's, so that we must only consider, say, m > when summing independent 
correlations. Nevertheless, even restraining to m > one double-counts two quantities, to wit 
and g^Q- Moreover, since g® is always zero, for a given I one can have a total of 6£ + 3 non-null 
independent quantities. Instead of restricting ourselves to m > (and have to worry about f ei j 
and gf^) we shall count these independent quantities in the following way: s <E {0, +1} and either 
m € {— £, £} (for s = 1) or m G {0, £} (for s = 0). With this convention the above relations simplifies 
to (using (B2)) 



(fern 

h 



Cm 



V2 



Ssl + Pz§sO 



Pv -M h 



V2 



c 



£+lm+s 



£m 3 

C t + (-l) s c 



s — s + ^1 



£+1 



(B15) 



and 



Jim J£' m' 
9im 9t> m' 



~ \ flm 
9lrn 



ft' m' / — fiei'Sss' {5mm' + SsoSmo) 



g}> 



-C e C e+1 6 u ,5 ss ,(5„ 



(B16) 



A similar estimator can be constructed in order to study the polarization. Following the same 
notation used in previous sections, we define 



f 



iXY _ 
im — 



1 



F, 



iXY 

£ m 



F, 



sXY 



1 

2~i 



zpsXY fjp. 
r £m ~ \ r £m 



;XY\ 
im J 



(B17) 
(B18) 



Applying then a Wigner rotation to Eq. (16), we find for such an estimator 

ft, 



(f£m Y ) 



<5 S 1 + ft^sO 



72 



sXY 
l tm J 



^ S UsXY 



(B19) 



and 



oxr = s-y c xy ,( ,y -s + X r XY 
n £m = C£+i m+s ^£ -r { *■) c £m W+l I 



rsXY rs'HK\ 
J£m J I'm' J 

sXY n s'HK\ 
k m i)£' m' J 





(ft 


HK 
m' 


~U X m Y ) 


u 


HK 
m' 



-Cf H Cj_^ 5u'5 ss > (5 rnm i + S s0 S mQ ) 
-Cf H Cj^f 5wS ss ' (S mm i — S S QS m o) 



(B20) 



where the coefficients c, 



-s±E,B 



are given by 

a ± E,B 



s±,T 



P 



(B21) 



and c s t ± T 

for I = 15 the difference is smaller than 1%. 



r, , = (see Eq. (B8)). Let us also note that for large £, c] 



-s±E,B 



c S p . In particular, 
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